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Comparison  of  Two  Methods  for  Estimating  Discharge 
and  Nutrient  Loads  from  Tidally  Affected  Reaches  of  the 
Myakka  and  Peace  Rivers,  West-Central  Florida 


By  Victor  A.  Levesque  and  K.M.  Hammett 


Abstract 

The  Myakka  and  Peace  River  Basins  consti¬ 
tute  more  than  60  percent  of  the  total  inflow  area 
and  contribute  more  than  half  the  total  tributary 
inflow  to  the  Charlotte  Harbor  estuarine  system. 
Water  discharge  and  nutrient  enrichment  have 
been  identified  as  significant  concerns  in  the  estu¬ 
ary,  and  consequently,  it  is  important  to  accurately 
estimate  the  magnitude  of  discharges  and  nutrient 
loads  transported  by  inflows  from  both  rivers. 

Two  methods  for  estimating  discharge  and 
nutrient  loads  from  tidally  affected  reaches  of  the 
Myakka  and  Peace  Rivers  were  compared.  The 
first  method  was  a  tidal-estimation  method,  in 
which  discharge  and  nutrient  loads  were  estimated 
based  on  stage,  water- velocity,  discharge,  and 
water-quality  data  collected  near  the  mouths  of  the 
rivers.  The  second  method  was  a  traditional  basin- 
ratio  method  in  which  discharge  and  nutrient  loads 
at  the  mouths  were  estimated  from  discharge  and 
loads  measured  at  upstream  stations. 

Stage  and  water-velocity  data  were  collected 
near  the  river  mouths  by  submersible  instruments, 
deployed  in  situ,  and  discharge  measurements 
were  made  with  an  acoustic  Doppler  current  pro¬ 
filer.  The  data  collected  near  the  mouths  of  the 
Myakka  River  and  Peace  River  were  filtered, 
using  a  low-pass  filter,  to  remove  daily  mixed-tide 
effects  with  periods  less  than  about  2  days.  The  fil¬ 
tered  data  from  near  the  river  mouths  were  used  to 
calculate  daily  mean  discharge  and  nutrient  loads. 
These  tidal-estimation-method  values  were  then 


compared  to  the  basin-ratio-method  values.  Four 
separate  30-day  periods  of  differing  streamflow 
conditions  were  chosen  for  monitoring  and  com¬ 
parison. 

Discharge  and  nutrient  load  estimates  com¬ 
puted  from  the  tidal-estimation  and  basin-ratio 
methods  were  most  similar  during  high-flow  peri¬ 
ods.  However,  during  high  flow,  the  values  com¬ 
puted  from  the  tidal-estimation  method  for  the 
Myakka  and  Peace  Rivers  were  consistently  lower 
than  the  values  computed  from  the  basin-ratio 
method.  There  were  substantial  differences 
between  discharges  and  nutrient  loads  computed 
from  the  tidal-estimation  and  basin-ratio  methods 
during  low-flow  periods.  Furthermore,  the  differ¬ 
ences  between  the  methods  were  not  consistent. 
Discharges  and  nutrient  loads  computed  from  the 
tidal-estimation  method  for  the  Myakka  River 
were  higher  than  those  computed  from  the  basin- 
ratio  method,  whereas  discharges  and  nutrients 
loads  computed  by  the  tidal-estimation  method  for 
the  Peace  River  were  not  only  lower  than  those 
computed  from  the  basin-ratio  method,  but  they 
actually  reflected  a  negative,  or  upstream,  net 
movement.  Short-term  tidal  measurement  results 
should  be  used  with  caution,  because  antecedent 
conditions  can  influence  the  discharge  and  nutrient 
loads.  Continuous  tidal  data  collected  over  a  1-  or 
2-year  period  would  be  necessary  to  more  accu¬ 
rately  estimate  the  tidally  affected  discharge  and 
nutrient  loads  for  the  Myakka  and  Peace  River 
Basins. 


Abstract  1 


INTRODUCTION 

The  Myakka  and  Peace  River  Basins  (fig.  1) 
constitute  more  than  60  percent  of  the  total  inflow  area 
and  contribute  more  than  half  the  total  tributary  inflow 
to  the  Charlotte  Harbor  estuarine  system  (Hammett, 
1990).  Nutrient  enrichment  has  been  identified  as  a 
significant  concern  in  the  estuary.  Consequently,  accu¬ 
rate  estimates  of  the  magnitude  of  nutrient  loads  trans¬ 
ported  to  Charlotte  Harbor  by  inflows  from  both  rivers 
are  important. 

Flows  from  approximately  62  percent  of  the 
Myakka  River  Basin  and  1 1  percent  of  the  Peace  River 
Basin  currently  (1996)  are  not  monitored  by  recording 
streamflow  gages  because  of  the  complications  associ¬ 
ated  with  computing  continuous  streamflow  in  tidally 
affected  reaches  of  the  rivers.  River  discharge  and 
constituent  loads  from  tidally  affected  areas  tradition¬ 
ally  have  been  estimated  using  a  drainage-area  or 
basin  ratio  to  extrapolate  data  from  upstream  stream- 
flow  gages  to  tidally  affected  downstream  stations. 
Richards  (1989)  evaluated  the  basin-ratio  method  for 
estimating  constituent  loads  from  unmonitored  inflow 
areas  around  lakes  and  determined  that  the  basin-ratio 
method  overestimated  the  actual  loads  more  than  half 
the  time.  Stoker,  Levesque,  and  Woodham  (1996)  also 
reported  differences  between  basin-ratio  extrapolated 
loads  and  loads  computed  from  measurements  in  the 
tidal  reaches  of  the  Alafia  and  Hillsborough  Rivers  in 
west-central  Florida. 

Technological  advances  have  made  it  possible 
to  measure  discharge  in  tidally  affected  reaches  with 
reasonable  accuracy.  Field  measurement  of  discharge 
makes  it  possible  to  estimate  continuous  discharge  in 
the  tidal  reaches  and  provides  another  method  of  inde¬ 
pendently  estimating  constituent  loads  (tidal-estima¬ 
tion  method). 

Accurate  estimates  of  discharge  and  constituent 
loads  from  the  Myakka  and  Peace  Rivers  are  needed 
so  that  water-resource  managers  can  make  informed 
long-term  monitoring  and  preservation  plans  for  the 
river  basins  and  the  Charlotte  Harbor  estuarine  sys¬ 
tem.  The  U.S.  Geological  Survey  (USGS)  began  a 
study  in  1994  to  measure  discharge  and  determine 
nutrient  loads  near  the  mouths  of  the  Myakka  and 
Peace  Rivers.  The  study  was  done  in  cooperation  with 
the  Southwest  Florida  Water  Management  District 
(SWFWMD)  and  was  jointly  funded  by  both  agencies. 
The  tidal-estimation  method  was  used  to  compute  dis¬ 
charge  and  nutrient  loads  from  data  collected  near  the 


mouths  of  the  rivers.  Discharge  and  nutrient  loads 
computed  from  the  tidal-estimation  method  were  then 
compared  to  those  computed  using  the  basin-ratio 
method. 

Purpose  and  Scope 

The  puipose  of  this  report  is  to  describe  and 
compare  two  methods — the  tidal-estimation  and 
basin-ratio  methods — that  were  used  for  estimating 
discharge  and  nutrient  loads  for  the  tidal  reaches  of  the 
Myakka  and  Peace  River  Basins.  Data  were  collected 
during  four  separate  30-day  periods  during  1994  and 
1995;  two  of  the  periods  represented  low-flow  condi¬ 
tions  (May-June  1994  and  April-May  1995)  and  two 
represented  high-flow  conditions  (September-October 
1994  and  October-November  1995).  Both  methods 
were  used  to  compute  discharge  and  nutrient  loads  for 
the  data-collection  periods,  and  the  results  were  then 
compared. 

Description  of  Study  Area 

The  Myakka  and  Peace  Rivers  empty  into 
northern  Charlotte  Harbor  (fig.  1).  The  Myakka  River 
receives  surface  runoff  from  an  area  of  approximately 
602  mi2.  The  headwaters  of  the  Myakka  River  are  in 
eastern  Manatee  County  near  Myakka  Head,  and  the 
river  flows  in  a  southerly  direction  through  Manatee, 
Sarasota,  and  Charlotte  Counties,  where  it  empties 
into  northwestern  Charlotte  Harbor.  Approximately 
373  mi2  (62  percent)  of  the  Myakka  River  Basin  drains 
into  the  river  downstream  from  the  continuous-record 
gaging  station  near  Sarasota  (station  5,  fig.  1). 

The  Peace  River  receives  surface  runoff  from  an 
area  of  approximately  2,350  mi2.  The  headwaters  of 
the  Peace  River  are  in  central  Polk  County,  and  the 
river  flows  in  a  southerly  direction  through  Polk, 
Hardee,  De  Soto,  and  Charlotte  Counties,  where  it 
empties  into  northeastern  Charlotte  Harbor.  Approxi¬ 
mately  260  mi2  (11  percent)  of  the  Peace  River  Basin 
drains  into  the  river  downstream  from  the  continuous- 
record  gaging  station  at  Arcadia  (station  1,  fig.  1)  and 
the  tributary  gages  on  Joshua  Creek  (station  2,  fig.  1), 
Horse  Creek  (station  3,  fig.  1),  and  Shell  Creek 
(station  4,  fig.  1). 
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Figure  1.  Myakka  River  and  Peace  River  Basins,  west-central  Florida. 
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METHODS 

Discharge-measurement  and  water-quality 
sampling  methods  used  for  this  study  were  adapted 
from  previously  established  techniques  (Simpson  and 
Oltmann,  1992;  Stoker,  Levesque,  and  Fritz,  1996). 
Reconnaissance  of  the  study  area  was  conducted  in 
February  1994  to  locate  tidal  stations  near  the  mouths 
of  the  Myakka  and  Peace  Rivers  that  would  be  suitable 
for  instrument  deployment  and  for  making  discharge 
measurements.  Bathymetric  and  discharge  characteris¬ 
tics  were  evaluated  to  determine  placement  of  record¬ 
ing  submersible  instrument  packages  (SIPs).  The  SIPs 
were  used  to  measure  and  record  instantaneous  water 
velocity,  water  temperature,  specific  conductance,  and 
water  depth  at  15 -minute  intervals  during  each  of  the 
30-day  deployment  periods.  The  four  deployments 
periods  were  May  3  to  June  2,  1994,  September  7  to 
October  6,  1994,  April  17  to  May  18,  1995,  and 
October  23  to  November  22,  1995.  During  the  deploy¬ 
ments,  tidal  station  discharge  measurements  were 
made  at  selected  times  near  the  mouths  of  the  rivers 
using  an  acoustic  Doppler  current  profiler  (ADCP). 

At  the  continuous-record-nontidal  stream  flow  stations 
upstream,  discharge  measurements  were  made  through¬ 
out  1994-95  using  standard  USGS  methods  (Rantz  and 
others,  1982a;  1982b).  Water-quality  samples  were 
collected  at  the  discharge-measurement  stations  near 
the  mouths  of  the  rivers  and  at  two  of  the  upstream 
gaging  stations  (stations  1  and  5,  fig.  1)  during  each 
deployment. 

Discharge 

Water-surface-elevation  (hereafter  referred  to  as 
stage),  water-velocity,  and  discharge  data  were  used  to 
establish  a  stage- velocity-discharge  relation  for  each 
tidal  station  (stations  6  and  7,  fig.  1).  A  relation  between 
stage  and  cross-sectional  area  was  determined  by 
using  a  fathometer  to  measure  depths  of  the  cross 
sections  where  discharge  measurements  were  made. 
Regression  analysis  was  used  to  relate  mean  velocity 
in  the  cross  section  (measured  with  the  ADCP)  to  an 
index  velocity  (measured  with  the  SIPs).  Discharge 
was  then  computed  as  the  product  of  cross-sectional 
area  and  mean  velocity  (Rantz  and  others, 
1982a;1982b). 

The  configuration  of  the  SIPs  is  described  in 
detail  by  Levesque  and  Schoellhamer  (1995).  Each 
SIP  consisted  of  an  aluminum  stand  with  two  electro¬ 
magnetic  current  meters  mounted  about  4  ft  above  the 


base  of  the  stand,  water  temperature  and  conductance 
sensors  mounted  about  3  ft  above  the  base,  and  a  pres¬ 
sure  transducer  mounted  near  the  base.  The  current 
meters,  conductance  sensor,  and  transducer  were  con¬ 
nected  to  an  electronic  data  logger  that  recorded  mea¬ 
surements.  Each  SIP  had  a  communication  cable  that 
allowed  the  data  logger  to  be  queried  and  checked 
from  the  surface. 

A  description  of  the  ADCP  is  provided  by  Sim¬ 
pson  and  Oltmann  (1992).  The  ADCP  boat-mounted 
instrument  acoustically  measures  water  velocities,  and 
accompanying  software  calculates  discharge  as  the 
boat  traverses  the  measurement  section.  Discharge 
measurements  in  this  study  were  made  over  a  range  of 
tidal  conditions  during  the  SIP  deployments. 

After  each  deployment,  SIP  data,  fathometer 
data,  and  ADCP  discharge  data  were  used  to  develop 
stage/area  and  index-velocity/mean-velocity  relations. 
All  stage  data  were  referenced  to  sea  level.  Instanta¬ 
neous  discharge  was  computed  for  15 -minute  intervals 
when  velocity  and  stage  data  were  available  during  the 
30-day  deployments. 

Discharges  in  tidal  rivers  have  a  high-frequency 
variation  that  dominates  the  discharge  patterns.  This 
variation  is  characterized  by  frequent  reversals  in  flow 
and  changes  in  flow  magnitude  and  is  caused  by  the 
upstream  and  downstream  movements  of  water  that 
result  from  tidal  currents.  The  amount  of  discharge 
caused  by  basin  inflow  can  be  small  relative  to  the 
amount  of  discharge  caused  by  tidal  currents.  However, 
the  basin  inflow  and  accompanying  nutrient  loads  are 
the  factors  of  interest  in  this  and  many  other  studies. 
Astronomical  tidal  variations  in  time-series  data,  such 
as  tidal  discharge,  can  be  removed  using  a  mathemati¬ 
cal  filter.  The  filtered  data  can  then  be  examined  to 
determine  variations  in  the  data  that  are  the  result  of 
longer  term  (greater  than  2  days)  astronomical  tidal 
variations  as  well  as  variations  resulting  from  nontidal 
processes,  such  as  basin  inflow  and  meteorological 
patterns.  The  Godin  filter,  a  low-pass  digital  filter 
(Walters  and  Heston,  1982),  was  used  to  remove  short¬ 
term  (less  than  about  2  days)  variations  in  instanta¬ 
neous  discharge  caused  by  diurnal  and  semidiurnal 
tides  (mixed  tides).  Once  mixed-tide  variations  were 
removed,  daily  mean  discharges  were  computed 
(Stoker,  Levesque,  and  Woodham,  1996).  Digitally 
filtering  the  discharge  data  results  in  a  loss  of  about 
2  days  of  data  near  the  beginning  and  end  of  a  data- 
collection  period.  Therefore,  a  30-day  period  of  data 
collection  results  in  digitally  filtered  data  for  a  period 
of  about  25  days. 


4  Comparison  of  Two  Methods  for  Estimating  Discharge  and  Nutrient  Loads  from  Tidally  Affected  Reaches  of  the  Myakka  and 
Peace  Rivers,  West-Central  Florida 


Standard  USGS  discharge-measurement  tech¬ 
niques  were  used  at  the  continuous-record  gaging 
stations  in  the  nontidal  reaches  of  both  rivers  during 
1994-95.  The  stations  measured  were  Myakka  River 
near-  Sarasota  (station  5,  fig.  1),  Peace  River  at  Arcadia 
(station  1,  fig.  1),  Joshua  Creek  at  Nocatee  (station  2, 
fig.  1),  Horse  Creek  near  Arcadia  (station  3,  fig.  1), 
and  Shell  Creek  near  Punta  Gorda  (station  4,  fig.  1). 


Myakka  River 

The  SIP,  near  the  mouth  of  the  Myakka  River 
(fig.  2),  was  positioned  about  1,500  ft  upstream  from 
the  El  Jobean  bridge  (State  Road  776)  and  about  1,000  ft 
from  the  south  bank.  A  few  days  after  the  SIP  was 
deployed,  ADCP  discharge  measurements  began,  and 
fathometer  measurements  were  made  at  the  measure¬ 
ment  section  downstream  from  the  El  Jobean  Bridge. 


Figure  2.  Location  of  submersible  instrument  package  and  discharge-measurement  section  near  the  mouth  of  the 
Myakka  River,  west-central  Florida. 
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The  measurement  cross  section  had  depths  ranging 
from  6  to  17  ft  and  an  average  depth  of  about  10  ft. 
The  section  was  about  1 ,200  ft  in  length  and  was  con¬ 
fined  on  both  edges  by  concrete  seawalls.  Individual 
discharge  measurements  required  about  10  minutes  to 
complete,  and  more  than  200  individual  measurements 
were  made  during  the  four  deployment  periods.  Three 
days  of  SIP  data  were  lost  during  the  September- 
October  1994  deployment  because  water  leaked  into 
the  communication  cable. 


Peace  River 

The  SIP  near  the  mouth  of  the  Peace  River  (fig.  3) 
was  positioned  approximately  equidistant  between  the 
north-  and  southbound  bridges  of  U.S.  Highway  41 
and  about  1,000  ft  from  the  south  bank.  ADCP  dis¬ 
charge  measurements  and  fathometer  transects  also 
were  made  at  a  cross  section  between  the  bridges. 

The  measurement  cross  section  was  characterized  by  a 
gradually  sloping  bottom  with  depths  ranging  from 
1  to  12  ft  and  an  average  depth  of  about  7  ft.  The  cross 


Figure  3.  Location  of  submersible  instrument  package  and  discharge-measurement  section  near  the  mouth  of  the 
Peace  River,  west-central  Florida. 


6  Comparison  of  Two  Methods  for  Estimating  Discharge  and  Nutrient  Loads  from  Tidally  Affected  Reaches  of  the  Myakka  and 
Peace  Rivers,  West-Central  Florida 


section  was  about  4,100  ft  in  length  with  concrete  sea¬ 
walls  at  both  ends,  and  approximately  700  ft  of  the 
measurement  section  was  not  measurable  due  to 
ADCP  limitations  in  water  depths  of  less  than  6  ft. 

The  700  ft  not  directly  measured  was  comprised  of 
two  sections,  approximately  500  ft  near  the  north  bank 
and  approximately  200  ft  near  the  south  bank.  The  dis¬ 
charges  were  estimated  for  the  sections  that  were  not 
measured  using  the  measured  discharge  and  the  ADCP 
software.  Discharges  in  the  sections  that  were  not 
measured  ranged  from  less  than  5  percent  to  about  10 
percent  of  the  total  discharge.  Individual  discharge 
measurements  required  about  20  minutes  to  complete, 
and  more  than  200  individual  measurements  were 
made  during  the  four  deployment  periods.  A  loss  of 
SIP  data  occurred  during  two  of  the  four  deployments. 
All  data  for  6  days  were  lost  during  the  May-June 
1994  deployment  as  a  result  of  water  leaking  into  the 
communications  cable,  and  9  days  of  velocity  data 
were  lost  during  the  October-November  1995  deploy¬ 
ment  as  a  result  of  water  leaking  into  the  SIP 

Water  Quality 

Water-quality  samples  were  collected  during 
ebb-,  slack-,  and  floodtides  at  the  ADCP  measurement 
sections  during  SIP  deployments.  Depth-integrated 
water-quality  samples  were  collected  using  a  weighted 
glass  bottle  and  were  transferred  into  a  high-density 
polyethylene  water-sample  churn  for  splitting  into 
individual  bottles  for  analyses.  Water-quality  samples 
were  collected  during  ADCP  measurements  on  the 
basis  of  flow  distributions  so  that  the  location  of  the 
sampling  was  in  the  center  of  sections  that  had  approx¬ 
imately  equal  discharge;  samples  were  collected  at 
three  stations  along  the  measurement  section  on  the 
tidal  Myakka  River  and  at  four  stations  along  the  mea¬ 
surement  section  on  the  tidal  Peace  River.  During 
water-quality  sampling,  field  measurements  of  specific 
conductance,  pH,  water  temperature,  and  dissolved 
oxygen  were  made  at  2-ft  intervals  to  describe  vertical 
variability  at  each  sampling  station.  Quality-assurance 
samples  (duplicates  or  blanks)  were  collected  and  sent 
with  the  regular  samples.  The  water-quality  data  col¬ 
lected  for  this  study  can  only  be  considered  represen¬ 
tative  of  the  four  deployment  periods  between  1994 
and  1995.  Long-term  data  would  be  required  to  accu¬ 
rately  characterize  nutrient  concentration  ranges  for 
the  tidal  stations  during  other  time  periods. 

Water-quality  samples  also  were  collected  at  the 
nontidal-recording  streamflow  stations  on  the  Myakka 
and  Peace  Rivers  (stations  5  and  1)  during  the  deploy¬ 


ments  to  enhance  the  routine  bimonthly  sampling  that 
is  done  during  the  year.  Depth-integrated  water-quality 
samples  were  collected  at  equally  spaced  distances  in 
the  cross  section.  Duplicate  or  blank  field  quality- 
assurance  samples  were  collected  and  sent  with  the 
regular  samples. 

Water  samples  were  packed  in  a  plastic  cooler 
with  ice  and  sent  to  the  USGS  laboratory  in  Ocala, 
Fla.,  for  analyses.  Samples  were  analyzed  for  unfil¬ 
tered  nutrients  (ammonia  nitrogen,  ammonia  plus 
organic  nitrogen,  nitrite  nitrogen,  nitrite  plus  nitrate 
nitrogen,  total  phosphorus,  and  orthophosphorus), 
total  suspended  solids,  turbidity,  color,  and  specific 
conductance.  The  analytical  methods  are  documented 
by  Fishman  and  Friedman  (1989). 

Computation  of  Loads 

Nutrient  loads  are  a  measure  of  the  amount  of 
total  nitrogen  and  total  phosphorus  that  is  transported 
by  water  during  a  specified  time  interval.  The  load  of  a 
constituent  is  expressed  as  a  weight  per  unit  time,  such 
as  pounds  per  day  or  tons  per  month.  The  nutrient  load 
at  any  given  time  can  be  determined  if  the  nutrient 
concentration  and  the  discharge  at  the  time  of  sam¬ 
pling  are  known.  In  this  report,  loads  computed  for 
each  sample  are  considered  measured  loads.  Nutrient 
loads  for  the  tidal  and  nontidal  stations  were  computed 
using  the  following  equation: 

L  =  C  x  Q  x  0.002697, 

where 

L  is  the  nutrient  load,  in  short  tons  per  day; 

C  is  the  concentration  of  total  nitrogen  or  total 
phosphorus,  in  milligrams  per  liter; 

Q  is  discharge,  in  cubic  feet  per  second;  and 
0.002697  is  a  conversion  factor. 

Simple-linear  and  multiple  regression  analyses 
of  instantaneous  discharge  and  measured  nitrogen  and 
phosphorus  loads  were  used  to  estimate  continuous 
nutrient  loads  at  the  tidal  stations  on  the  Myakka  and 
Peace  Rivers  and  at  the  nontidal-upstream  stations. 
Because  load  is  zero  when  discharge  is  zero,  the 
regression  intercepts  were  adjusted  to  zero  (Helsel  and 
Hirsch,  1992).  Nutrient  loads  could  not  be  estimated 
for  periods  when  discharge  data  were  not  available. 

Tidal  Stations 

Nutrient  loads  measured  at  the  tidal  stations 
during  the  four  deployment  periods  were  used  to 
develop  regression  equations  that  described  the 
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relation  between  instantaneous  discharge  and  total 
nitrogen  and  total  phosphorus  load.  The  regression 
equations  were  then  applied  to  instantaneous  dis¬ 
charge  to  calculate  instantaneous  total  nitrogen  and 
total  phosphorus  loads  for  the  four  deployments. 

The  regression  equations  for  total  nitrogen  and  total 
phosphorus  loads  for  the  Myakka  River  at  El  Jobean 
(station  6.  fig.  1)  and  the  Peace  River  at  Punta  Gorda 
(station  7,  fig.  1)  are  shown  in  table  1. 


The  instantaneous  nitrogen  and  phosphorus 
loads  were  filtered  using  the  Godin  filter  to  remove 
short-term  mixed-tide  valuations  of  less  than  about 
2  days.  The  filtered  nutrient  loads  were  averaged  for 
daily  (24-hour)  periods,  and  these  daily  mean  values 
were  summed  to  compute  nutrient  loads  for  the  peri¬ 
ods  when  discharge  data  were  collected  at  the  tidal 
stations. 

Nontidal  Stations 

Prior  to  calculating  loads  for  the  upstream 
stations,  the  long-term  nutrient  concentration  and 
discharge  data  were  analyzed  for  trends  utilizing  the 
nonparametric  Kendall  Tau  test  (Smith  and  others, 
1982).  Trends  were  not  detected  in  total  nitrogen  and 
total  phosphorus  concentrations  (significance  level  of 
0.05)  for  any  station  during  the  periods  used  to 
develop  regression  equations.  Trends  were  not 
detected  in  annual  or  seasonal  (quarterly)  discharge 
for  the  Myakka  River  (station  5,  fig.  1  ),  Peace  River 
(station  1,  fig.  1  )  or  Horse  Creek  (station  3,  fig.  1  ). 
Nutrient  loads  for  the  Myakka  River,  Peace  River,  and 
Horse  Creek  stations  were  based  on  discharge  and 
water-quality  data  from  1986  to  1996.  Nutrient  loads 
for  Joshua  Creek  (station  2,  fig.  1)  were  based  on  dis¬ 
charge  and  water-quality  data  from  1991  to  1996,  and 
loads  for  Shell  Creek  (station  4,  fig.  1)  were  based  on 


discharge  and  water-quality  data  from  1994  to  1995. 
The  regression  equations  used  to  calculate  nutrient 
loads  for  the  nontidal  stations  are  shown  in  table  2. 

Estimates  of  discharge,  total  nitrogen  load,  and 
total  phosphorus  load  from  ungaged  parts  of  a  river 
basin  are  often  made  by  extrapolating  loads  from  the 
gaged  part  of  the  basin.  For  example,  if  50  mi2  of  a 
100-mi2  river  basin  were  gaged,  and  if  water-quality 
data  were  available  at  the  gaged  station  to  calculate 

nutrient  loads,  then  the 
total  river-basin  dis¬ 
charge  and  nutrient  loads 
could  be  estimated  by 
multiplying  discharge 
and  nutrient  loads  at  the 
gaged  station  by  2.00. 
The  basin-ratio  method 
assumes  that  the  basin 
yield  for  discharge  and 
nutrients  is  the  same  for 
the  ungaged  basin  as  the 
gaged  part  of  the  basin. 
Runoff,  ground-water  inflow,  and  nutrient  concentra¬ 
tions  can  vary  throughout  a  river  basin,  causing  an 
error  in  estimations  based  on  the  basin-ratio  method. 
For  this  study,  the  basin-ratio  method  was  applied  by 
using  the  daily  means  of  discharge,  total  nitrogen  load, 
and  total  phosphorus  load  for  the  nontidal  gaged  sta¬ 
tions.  The  total  nontidal  gaged  values  for  discharge, 
total  nitrogen  load,  and  total  phosphorus  load  were 
then  multiplied  by  a  ratio  to  account  for  the  ungaged 
part  of  each  basin.  Using  the  discharge  and  nutrient 
concentrations  from  station  5  for  the  Myakka  River 
and  stations  1,  2,  3,  and  sometimes  station  4  for  the 
Peace  River,  estimates  of  total  basin  discharge  and 
loading  were  computed  by: 

Basin-ratio  value  at  Myakka  River  at  El  Jobean  = 
Station  5  value  x  2.63 

For  deployments  in  1994, 

Basin-ratio  value  at  Peace  River  at  Punta  Gorda  = 
(Station  1  value  +  Station  2  value  +  Station  3  value)  x  1 .37 

For  deployments  in  1995  when  data  for  Shell  Creek 
were  available, 

Basin-ratio  value  at  Peace  River  at  Punta  Gorda  = 
(Station  1  value  +  Station  2  value  +  Station  3  value  + 
Station  4  value)  x  1.12. 


Table  1.  Regression  equations  used  for  tidal-estimation  nutrient  loads 

[Station  locations  are  shown  in  fig.  1.  R2,  coefficient  of  determination;  n,  number  of  observations;  Q, 
instantaneous  discharge,  in  cubic  feet  per  second] 


Station 

Dependent  variable 
(short  tons  per  day) 

R2 

Standard 

error 

(short  tons) 

n 

Equation 

Myakka  River  at  El  Jobean 

nitrogen  load  (N) 

0.98 

5.16 

20 

N  =  2.693xl0-3  x  Q 

(station  6) 

phosphorus  load  (P) 

.94 

2.08 

19 

P  =  6.547x10 4  x  Q 

Peace  River  at  Punta  Gorda 

nitrogen  load  (N) 

.98 

9.84 

16 

N  =  3.160xl0'3  x  Q 

(station  7) 

phosphorus  load  (P) 

.97 

4.59 

16 

P=  l.lllxlCHxQ 
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Table  2.  Regression  equations  used  for  basin-ratio  nutrient  loads 

[Station  locations  are  shown  in  fig.  1 .  R2,  coefficient  of  determination;  n,  number  of  observations;  Q,  instantaneous  discharge,  in 
cubic  feet  per  second] 


Station 

Dependent  variable 
(short  tons  per  day) 

R2 

Standard 

error 

(short  tons) 

n 

Equation 

Myakka  River  near  Sarasota 

nitrogen  load  (N) 

0.98 

0.25 

32 

N  =  (3.628x10 3)  x  Q  -  (4.000  xlO'7)  x  Q2 

(station  5) 

phosphorus  load  (P) 

.98 

.13 

32 

P  =  (1.437x10  3)xQ 

Peace  River  at  Arcadia 

nitrogen  load  (N) 

.98 

.88 

63 

N  =  (5.324X10'3)  x  Q  -  (2.000  x  10 7)  x  Q2 

(station  1) 

phosphorus  load  (P) 

.97 

.51 

63 

P  =  (2.671x10 3)  x  Q  -  (9.520  x  10  *)  x  Q2 

Joshua  Creek  at  Nocatee 

nitrogen  load  (N) 

.93 

.21 

21 

N  =  (6.164xlO'3)xQ 

(station  2) 

phosphorus  load  (P) 

.98 

.02 

21 

P  =  (7.623x10 4)  x  Q 

Horse  Creek  near  Arcadia 

nitrogen  load  (N) 

.98 

.26 

57 

N  =  (4.585xl0'3)  x  Q  -  (4.000x10 7)  x  Q2 

(station  3) 

phosphorus  load  (P) 

.98 

.11 

57 

P  =  (1.495x10 3)  x  Q  -  (6.344x10  s)  x  Q2 

Shell  Creek  near  Punta  Gorda 

nitrogen  load  (N) 

.99 

.27 

7 

N  =  (2.965X10'3)  x  Q 

(station  4) 

phosphorus  load  (P) 

.99 

.06 

7 

P  =  (4.980x10 4)  x  Q 

COMPARISON  OF  METHODS  USED  TO 
ESTIMATE  DISCHARGE  AND  NUTRIENT 
LOADS 

Discharges  and  nutrient  loads  for  the  Myakka 
and  Peace  River  Basins  were  estimated  by  the  tidal- 
estimation  method  and  the  basin-ratio  method.  Several 
factors  affect  the  relative  accuracy  of  each  of  the  meth¬ 
ods.  The  basin-ratio  method  assumes  that  the  quality, 
quantity,  and  time  distribution  of  runoff  and  base  flow 
in  the  tidally  affected  reaches  are  the  same  as  they  are 
in  the  upstream  paid  of  the  basin.  However,  the  upland 
parts  of  the  Myakka  and  Peace  River  Basins  have 
steeper  land-surface  slopes  than  the  coastal  areas;  urban 
development  tends  to  be  centered  in  the  ungaged 
coastal  areas,  whereas  agricultural  development  pre¬ 
dominates  in  the  upland  areas.  Thus,  many  of  the  con¬ 
ditions  that  affect  the  quality,  quantity,  and  time 
distribution  of  discharge  are  different  in  the  tidally 
affected  reaches  than  in  the  upland  reaches.  The  basin- 
ratio  method  also  assumes  that  discharge  is  hydrauli¬ 
cally  controlled  at  the  upstream  end  of  the  reach,  and 
does  not  take  into  account  any  effects  that  tides  might 
have  on  the  hydraulic  characteristics  of  flow.  In  west- 
central  Florida,  where  stream-channel  gradients  are 
low,  velocities  associated  with  stream  discharge  are 
also  low  and  tidal  forces  have  a  major  impact  on  flow 
characteristics  near  the  mouths  of  coastal  rivers. 

ADCP  technology  has  greatly  enhanced  the  ability  of 
hydrologists  to  make  field  measurements  of  instanta¬ 


neous  discharge  in  estuarine  areas.  However,  the  tidal- 
estimation  method  is  also  subject  to  errors  that  can  be 
significant  relative  to  the  magnitude  of  daily  mean 
dischai'gc. 

Myakka  River 

More  than  200  discharge  measurements  were 
made  using  the  ADCP  at  the  Myakka  River  at  El  Jobean 
(fig.  2).  Field  measurements  did  not  reveal  any  unusual 
flow  characteristics  in  the  vicinity  of  the  measurement 
section  or  near  the  location  of  the  submersible  instru¬ 
ment  package.  Ebbtide  and  floodtide  discharges  were 
relatively  uniform  across  the  measurement  section. 

At  slacktide,  the  dischai'gc  typically  had  bidirectional 
characteristics,  with  water  moving  upstream  at  one  side 
of  the  measurement  section  and  downstream  at  the 
other  side.  Unidirectional  discharges  recurred  within 
about  20  minutes  following  slacktide. 

Instantaneous  discharge  values  were  computed 
using  the  index/mean  velocity  rating  and  the 
stage/area  rating  and  ranged  from  about  -20,000  to 
+20,000  ft3/s  during  the  four  deployment  periods. 
Instantaneous  discharges  for  May  1994  arc  shown  in 
figure  4a.  In  this  report,  negative  discharge  values 
denote  discharge  moving  in  an  upstream  direction. 

The  standard  error  in  instantaneous  discharge,  which 
was  computed  by  multiplying  the  average  standard 
error  of  the  index- velocity /mean- velocity  regressions 
by  the  average  cross-sectional  area  of  the  measure¬ 
ment  section,  averaged  about  ±  1,100  ft3/s. 
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DISCHARGE,  DISCHARGE,  DAILY  MEAN  BASIN  DISCHARGE, 

IN  CUBIC  FEET  PER  SECOND  IN  CUBIC  FEET  PER  SECOND  RAINFALL,  IN  INCHES  IN  CUBIC  FEET  PER  SECOND 
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Figure  4.  Instantaneous  discharge,  tidal-estimation  discharge,  basin-ratio  discharge,  daily  mean  basin 
rainfall,  and  stage  for  the  Myakka  River  at  El  Jobean,  west-central  Florida,  May  3  to  June  1 ,  1 994. 
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STAGE, 

IN  FEET  ABOVE  SEA  LEVEL 


The  Godin  filter  was  used  to  remove  short-term 
variations  in  instantaneous  discharge  and  stage  that 
were  caused  by  tides  (fig.  4a,  c,  d).  Removal  of  the 
short-term  valuations  makes  it  possible  to  examine 
variations  that  occur  over  periods  longer  than  2  days. 
Filtered  discharge  data  typically  is  out  of  phase  with 
filtered  stage  data  (fig.  4d),  indicating  that  when  stage 
is  high,  as  it  was  from  May  5-10,  discharge  is  sup¬ 
pressed.  Conversely,  when  the  stage  declines  over 
several  days,  as  it  did  from  May  15-20,  discharge 
increases.  Stage,  therefore,  has  a  significant  effect  on 
dischai'gc  during  particular  conditions.  After  varia¬ 
tions  resulting  from  short-term  tidal  action  were 
removed,  discharge  ranged  from  about  - 1 ,000  ft  Vs  in 
November  1995  to  about  5,000  ft3/s  in  September 
1994.  Daily  mean  discharge,  computed  from  the  fil¬ 
tered  dischai'gc  data  for  Myakka  River  at  El  Jobean 
(fig.  2),  ranged  from  an  average  of  649  ft3/s  during 
April-May  1995,  to  an  average  of  4,079  ft3/s  during 
September-October  1994. 

Daily  mean  values  of  basin  rainfall,  discharge, 
specific  conductance,  total  nitrogen  load,  and  total 
phosphorus  load  for  the  Myakka  River  at  El  Jobean 
are  shown  in  figures  5  through  8  for  the  four  deploy¬ 
ment  periods.  Rainfall  for  the  basin  was  estimated  by 
averaging  rainfall  at  the  National  Weather  Service 
stations  at  Myakka  River  State  Park  and  Fort  Green 
(National  Oceanic  and  Atmospheric  Administration, 
1994-1995).  The  rainfall  data  do  not  appeal-  to  be  a 
good  indicator  of  discharge  at  the  mouth  of  the 
Myakka  River  and  are  probably  the  result  of  the 
number  and  location  of  the  rainfall  stations  which 
were  used  to  compute  daily  rainfall  for  the  basin. 

Data  shown  in  figures  5-8  illustrate  the  effect 
that  downstream  stage  has  on  controlling  discharge 
near  the  mouth  of  the  Myakka  River  over  all  four 
deployments.  Stage  and  discharge  are  generally 
inversely  related  near  the  mouth  of  the  river.  As  stage 
decreases,  an  increase  in  discharge  follows,  and  as 
stage  increases,  discharge  decreases. 

Generally,  discharge  and  specific  conductance 
also  are  inversely  related.  The  lowest  specific  conduc¬ 
tance  occurred  during  September-October  1994 
(fig.  6c)  when  discharges  were  the  highest.  The  high¬ 
est  specific  conductance  values  (greater  than  36,000 
pS/cm)  were  measured  during  May-June  1994 
(fig.  5c).  Discharges  were  actually  lower  during  the 
April-May  1995  (fig.  7b)  than  during  May-June  1994, 
but  stage  near  the  mouth  was  higher  during  the  1994 
deployment  and,  combined  with  low  discharges, 


resulted  in  the  highest  specific  conductance  values. 
The  higher  stages  during  the  May-June  1994  deploy¬ 
ment  represented  a  greater  upstream  movement  of 
saline  water  from  Charlotte  Harbor,  as  indicated  by  the 
higher  specific  conductance  values  (fig.  5c). 

Using  discharges  computed  from  the  tidal- 
estimation  method,  daily  mean  nutrient  loads  at  the 
Myakka  River  at  El  Jobean  were  negative  for  some 
days  during  October-November  1995  (fig.  8d),  indi¬ 
cating  that  at  times,  nutrients  were  going  into  storage 
in  the  reach  upstream  from  the  mouth.  The  net  loads 
for  all  four  deployment  periods  were  positive,  indicat¬ 
ing  that  overall  there  was  a  transport  of  nutrients  into 
Charlotte  Harbor  during  each  deployment. 

Total  discharge  and  total  nutrient  loads  com¬ 
puted  for  each  deployment  are  presented  in  table  3. 
Using  the  tidal-estimation  method,  total  nitrogen  loads 
at  the  mouth  of  the  Myakka  River  ranged  from  44 
short  tons  during  October-November  1995  to  143 
short  tons  during  September-October  1994.  Using  the 
basin-ratio  method,  total  nitrogen  loads  ranged  from 
12  short  tons  for  April-May  1995  to  193  short  tons  for 
September-October  1994.  Using  the  tidal-estimation 
method,  total  phosphorus  loads  ranged  from  1 1  short 
tons  during  April-May  and  October-November  1995 
to  59  short  tons  during  September-October  1994.  Total 
phosphorus  loads  computed  using  the  basin  ratio 
method  ranged  from  5  short  tons  for  April-May  1995 
to  99  short  tons  for  September-October  1994. 

Discharges  and  nutrient  loads  computed  using 
the  tidal-estimation  method  for  the  two  low-flow 
periods  (May-June  1994  and  April-May  1995)  were 
greater  than  those  computed  using  the  basin-ratio 
method.  The  tidal-estimation  discharge  was  2.88  and 
5.22  times  greater  than  the  basin-ratio  discharge.  The 
tidal-estimation  nitrogen  load  was  2.18  and  3.89  times 
greater  than  the  basin-ratio  nitrogen  load,  and  the 
tidal-estimation  phosphorus  load  was  1.31  and  2.34 
times  greater  than  the  basin-ratio  phosphorus  load. 
Nutrient  concentrations  at  the  mouth  of  the  Myakka 
River  are  different  than  the  concentrations  upstream, 
which  is  evident  from  the  regression  equations  in 
tables  1  and  2.  If  nutrient  concentrations  remained 
unchanged  throughout  the  basin,  then  the  ratios  of 
tidal-estimation/basin-ratio  discharges  and  nutrient 
loads  would  be  the  same.  Because  nutrient  concentra¬ 
tions  change,  loads  computed  from  the  basin-ratio 
method  would  be  in  error,  even  if  discharges  com¬ 
puted  from  the  basin-ratio  method  were  completely 
accurate. 
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Figure  5.  Daily  mean  basin  rainfall,  stage,  tidal-estimation  discharge,  basin-ratio  discharge,  specific  conductance, 
tidal-estimation  nutrient  loads,  and  basin-ratio  nutrient  loads  for  the  Myakka  River  at  El  Jobean,  west-central 
Florida,  May  3  to  June  1 ,  1 994. 
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DAILY  MEAN  SPECIFIC  CONDUCTANCE,  DAILY  MEAN  STAGE,  IN  FEET 
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Figure  6.  Daily  mean  basin  rainfall,  stage,  tidal-estimation  discharge,  basin-ratio  discharge,  specific  conductance, 
tidal-estimation  nutrient  loads,  and  basin-ratio  nutrient  loads  for  the  Myakka  River  at  El  Jobean,  west-central  Florida, 
September  7  to  October  5,  1994. 
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Figure  7.  Daily  mean  basin  rainfall,  stage,  tidal-estimation  discharge,  basin-ratio  discharge,  specific  conductance, 
tidal-estimation  nutrient  loads,  and  basin-ratio  nutrient  loads  for  the  Myakka  River  at  El  Jobean,  west-central  Florida, 
April  17  to  May  17,  1995. 
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Figure  8.  Daily  mean  basin  rainfall,  stage,  tidal-estimation  discharge,  basin-ratio  discharge,  specific  conductance, 
tidal-estimation  nutrient  loads,  and  basin-ratio  nutrient  loads  for  the  Myakka  River  at  El  Jobean,  west-central  Florida, 
October  23  to  November  21,1 995. 
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Table  3.  Total  discharges  and  nutrient  loads  at  the  Myakka  River  at  El  Jobean,  computed  using 
the  tidal-estimation  and  basin-ratio  methods 

[Mgal  =  million  gallons;  ft3/s  =  cubic  feet  per  second] 


Deployment  period 

Element 

Tidal 

Basin-ratio 

Tidal 

Corrected 

basin 

ratio2 

estimation 

estimation1 

Basin  ratio 

May  3-  June  1,  1994 

Volume  (Mgal) 

13,954 

4,844 

2.88 

7.58 

(25  days  of  usable  data) 

Discharge  (ft3/ s) 

864 

300 

Nitrogen  Load  (short  tons) 

58 

27 

2.18 

5.73 

Phosphorus  Load  (short  tons) 

14 

11 

1.31 

3.46 

Sept.  7  -  Oct.  5,  1994 

Volume  (Mgal) 

34,269 

44,328 

.77 

2.03 

(13  days  of  usable  data) 

Discharge  (ft3/s) 

4,079 

5.276 

Nitrogen  Load  (short  tons) 

143 

193 

.74 

1.94 

Phosphorus  Load  (short  tons) 

59 

99 

.60 

1.57 

April  17  -  May  17,  1995  Volume  (Mgal) 

10,903 

2.091 

5.22 

13.72 

(26  days  of  usable  data) 

Discharge  (ft3/s) 

649 

124 

Nitrogen  Load  (short  tons) 

45 

12 

3.89 

10.24 

Phosphorus  load  (short  tons) 

11 

5 

2.34 

6.15 

Oct.  23 -Nov.  21,  1995 

Volume  (Mgal) 

10,684 

14,386 

.74 

1.95 

(25  days  of  usable  data) 

Discharge  (ft3/ s) 

661 

890 

Nitrogen  Load  (short  tons) 

44 

78 

.57 

1.51 

Phosphorus  Load  (short  tons) 

11 

32 

.33 

.88 

1  Basin-ratio  estimation  =  value  at  Station  5  (fig.  1)  *  2.63. 

2  Corrected  basin  ratio  *  value  at  Station  5  (fig.  1 )  =  Tidal  estimation. 


Computations  for  the  two  high-flow  deploy¬ 
ments  produced  different  results  than  for  the  low-flow 
periods  and  indicated  the  effects  of  seasonal  as  well  as 
spatial  differences  in  nutrient  concentrations.  For  the 
two  high-flow  periods,  discharges  and  nutrient  loads 
computed  using  the  tidal-estimation  method  were  less 
than  those  computed  using  the  basin-ratio  method, 
which  is  the  opposite  of  the  what  occurred  during  the 
low-flow  deployments.  Tidal-estimation  discharges 
were  less  than  basin-ratio  discharges  by  factors  of  0.77 
and  0.74.  Tidal-estimation  nitrogen  loads  were  lower 
than  basin-ratio  loads  by  factors  of  0.74  and  0.57  and 
tidal-estimation  phosphorus  loads  were  lower  by  fac¬ 
tors  of  0.60  and  0.33. 

Differences  in  estimates  of  discharge  and  nutri¬ 
ent  loads  computed  from  the  tidal-estimation  and 
basin-ratio  methods  can  result  from  several  factors.  As 
noted  previously,  topography,  land  use,  and  nutrient 
concentrations  are  different  in  the  gaged  upland  areas 
than  in  the  ungaged  coastal  areas,  and  the  basin-ratio 
method  does  not  take  these  differences  into  account. 
The  average  standard  error  in  measured  instantaneous 
discharges,  ±  1,100  ft3/s,  is  significant  relative  to  the 
magnitude  of  daily  mean  discharge,  particularly  dur¬ 
ing  low  flow.  The  errors  associated  with  index/mean 


velocity  regressions  tend  to  increase  as  flow  velocities 
decrease.  Consequently,  the  errors  in  the  tidal-estima¬ 
tion  method  tend  to  be  greater  during  low-flow  periods. 
The  basin-ratio  method  also  tends  to  have  greater 
errors  during  low-flow  periods,  because  stage  near  the 
mouth  of  the  river  has  an  increased  hydraulic  effect  on 
discharge  as  velocities  decrease;  the  basin-ratio 
method  has  no  way  of  accounting  for  this  effect. 

Differences  in  ground-water  inflow  in  the 
upland  and  tidal  reaches  may  contribute  to  differences 
in  discharge  and  nutrient  loads  computed  from  the  two 
methods.  Based  on  data  collected  between  1983  and 
1987,  Flammett  (1992)  reported  a  potential  for  move¬ 
ment  of  water  from  the  surficial  aquifer  system  into 
the  Myakka  River  at  the  upstream  end  of  the  tidal 
reach,  but  also  noted  there  was  no  detectable  flow  in 
the  river  during  low-flow  periods.  Potentiometric- 
surface  maps  of  the  intermediate  and  Floridan  aquifer 
systems  for  May  1994,  September  1994,  May  1995, 
and  September  1995  also  show  a  potential  for  ground- 
water  discharge  into  the  tidal  reach  of  the  Myakka 
River  (Metz  and  Brendle,  1994;  Metz,  1995;  Metz  and 
Stelman,  1995;  Metz  and  others  1996).  On  April  17, 
1995,  water  was  observed  flowing  from  a  rusted  well 
pipe  into  the  Myakka  River  downstream  of  the 
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El  Jobean  bridge.  The  water  had  a  strong  hydrogen 
sulfide  smell  indicative  of  ground  water,  further  evi¬ 
dence  of  possible  ground-water  seepage  into  the  tidal 
reach  of  the  Myakka  River.  Potentiometric-surface 
data  indicate  a  greater  potential  for  ground- water 
inflow  in  the  tidal  reach  than  in  the  area  upstream  from 
the  long-term  gaging  station  (station  5,  fig.  1),  which 
could  account  for  some  of  the  differences  in  discharge 
and  nutrient  loads  computed  from  the  tidal-estimation 
and  basin-ratio  methods. 

Another  factor  that  contributes  to  differences 
between  the  tidal-estimation  and  basin-ratio  estimates 
of  discharge  and  nutrient  loads  is  the  effect  of  water 
going  into  or  being  released  from  storage  in  the  tidal 
reach.  The  surface  area  of  the  tidal  reach  of  the 
Myakka  River  upstream  from  the  bridge  at  El  Jobean 
is  about  6  mi2.  If  inflow  to  the  tidal  reach  exceeded 
outflow  by  about  20  ft3/s  over  the  course  of  1  day,  the 
stage  in  the  reach  would  rise  just  0.01  ft.  Inflow  would 
have  to  exceed  outflow  by  about  200  ft3/s  over  the 
course  of  1  day  in  order  for  the  stage  in  the  reach  to 
rise  0. 1  ft.  Typical  tidal  range  at  the  Myakka  River  at 
El  Jobean  is  2.20  ft  (Hammett,  1992),  and  conse¬ 
quently,  a  change  in  stage  of  0. 1  ft  would  not  be  partic¬ 
ularly  noteworthy.  Water  will  go  into  or  be  released 
from  storage  depending  on  the  tidal  conditions  at  the 
downstream  end  of  the  tidal  reach  and  the  volume  of 
freshwater  discharge.  When  upstream  discharges  are 
low,  the  storage  factor  will  be  most  significant  and  the 
magnitude  and  timing  of  discharges  at  the  mouth 
could  be  substantially  different  than  discharges  esti¬ 
mated  from  the  upstream  gaging  station  using  the 
basin-ratio  method.  Since  deployments  for  this  study 
were  for  30-day  periods  separated  by  several  months, 
it  is  possible  that  for  the  two  low-flow  deployments, 
discharges  measured  with  the  ADCP  and  computed 
using  the  tidal-estimation  method  reflect  water  coming 
out  of  storage  in  the  tidal  reach.  By  collecting  continu¬ 
ous  data  for  a  longer  period  of  time  it  would  be  possi¬ 
ble  to  more  accurately  quantify  storage  effects, 
identify  seasonal  storage  variations,  and  reduce  the 
number  of  variables  contributing  to  errors  associated 
with  the  methods. 

Peace  River 

More  than  200  discharge  measurements  were 
made  using  the  ADCP  at  the  Peace  River  at  Punta 
Gorda  (fig.  3).  Analysis  of  the  field  measurements  did 
not  reveal  any  unusual  flow  characteristics  in  the 


vicinity  of  the  measurement  section  or  upstream  of  the 
U.S.  41  bridges.  The  ebbtide  and  floodtide  discharges 
were  relatively  uniform  across  the  measurement  sec¬ 
tion.  The  discharges  were  typically  bidirectional  at 
slacktide,  but  returned  to  unidirectional  within  about 
30  minutes  following  slacktide. 

Instantaneous  discharge  at  the  Peace  River  at 
Punta  Gorda  ranged  from  about  -40,000  to  more  than 
+40,000  ft3/s  during  the  four  deployment  periods. 
Instantaneous  discharges  for  April-May  1995  are 
shown  in  figure  9a.  The  standard  error  in  instanta¬ 
neous  discharge  averaged  about  ±  3,600  ft3/s.  After 
short-term  variations  in  instantaneous  discharge  were 
removed  using  the  Godin  filter,  discharge  ranged  from 
about  -2,000  ft3/s  in  April  and  November  1995  to 
about  +19,000  ft3/s  in  September  1994.  The  effects  of 
stage  on  discharge  in  the  Peace  River  during  some 
periods,  (fig.  9d)  are  similar  to  effects  in  the  Myakka 
River.  For  example,  during  May  6-10,  1995,  as  stage 
increased  discharge  decreased.  At  other  times  there 
appears  to  be  a  time  lag  between  decreasing  stage  and 
increasing  discharge. 

Daily  mean  values  of  basin  rainfall,  stage,  dis¬ 
charge,  specific  conductance,  total  nitrogen  load,  and 
total  phosphorus  load  for  the  Peace  River  at  Punta 
Gorda  are  shown  in  figures  10  through  13  for  the  four 
deployment  periods.  Rainfall  for  the  basin  was  esti¬ 
mated  by  averaging  rainfall  from  the  National  Weather 
Service  stations  (fig.  1)  at  Bartow,  Wauchula,  Arcadia, 
and  Punta  Gorda  (National  Oceanic  and  Atmospheric 
Administration,  1994-1995).  The  rainfall,  as  com¬ 
puted,  does  not  appear  to  be  a  good  indicator  of  total 
discharge  at  the  mouth  of  the  Peace  River.  The  lack  of 
correlation  between  computed  rainfall  and  discharge 
probably  is  related  to  the  timing  and  distribution  of 
rainfall  throughout  the  basin  as  well  as  the  small  num¬ 
ber  of  stations  that  were  used  to  compute  daily  mean 
rainfall. 

Data  shown  in  figures  10-13  illustrate  the  effect 
that  downstream  stage  has  on  discharge  near  the 
mouth  of  the  Peace  River  over  all  four  deployments. 
Downstream  stage  has  a  significant  hydraulic  effect  on 
discharge,  although  the  relation  between  stage  and  dis¬ 
charge  appears  to  be  more  complex  than  at  the  mouth 
of  the  Myakka  River.  In  the  tidal  reach  of  the  Peace 
River,  the  relative  timing  of  inflow  from  three  major 
tributaries  and  the  main  channel  affect  the  shape  of  the 
hydrograph  at  the  mouth.  Also,  there  is  a  lag  between 
the  time  a  storm  event  is  measured  at  the  upstream 
stations  and  the  time  the  storm  discharge  reaches  the 
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Figure  9.  Instantaneous  discharge,  tidal-estimation  discharge,  basin-ratio  discharge,  daily  mean  basin  rainfall,  and 
stage  for  the  Peace  River  at  Punta  Gorda,  west-central  Florida,  April  17  to  May  17,  1995. 
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Figure  10.  Daily  mean  basin  rainfall,  stage,  tidal-estimation  discharge,  basin-ratio  discharge,  specific  conductance, 
tidal-estimation  nutrient  loads,  and  basin-ratio  nutrient  loads  for  the  Peace  River  at  Punta  Gorda,  west-central  Florida, 
May  3  to  June  1,  1994. 
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Figure  11.  Daily  mean  basin  rainfall,  stage,  tidal-estimation  discharge,  basin-ratio  discharge,  specific  conductance, 
tidal-estimation  nutrient  loads,  and  basin-ratio  nutrient  loads  for  the  Peace  River  at  Punta  Gorda,  west-central  Florida, 
September  7  to  October  5,  1994. 
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Figure  12.  Daily  mean  basin  rainfall,  stage,  tidal-estimation  discharge,  basin-ratio  discharge,  specific  conductance, 
tidal-estimation  nutrient  loads,  and  basin-ratio  nutrient  loads  for  the  Peace  River  at  Punta  Gorda,  west-central  Florida, 
April  17  to  May  17,  1995. 
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Figure  13.  Daily  mean  basin  rainfall,  stage,  tidal-estimation  discharge,  basin-ratio  discharge,  specific  conductance,  tidal- 
estimation  nutrient  loads,  and  basin-ratio  nutrient  loads  for  the  Peace  River  at  Punta  Gorda,  west-central  Florida,  October  23 
to  November  21,  1995. 
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Table  4.  Total  discharges  and  nutrient  loads  at  the  Peace  River  at  Punta  Gorda,  computed  using  the 
tidal-estimation  and  basin-ratio  methods 

[Mgal  =  million  gallons;  ft3/s  =  cubic  feet  per  second] 


Deployment  period 

Element 

Tidal 

estimation 

Basin-ratio 

estimation1 

Tidal 

Corrected 

basin 

ratio2 

Basin  ratio 

May  3- June  1,  1994 

Volume  (Mgal) 

-2,396 

809 

-2.96 

-4.06 

(13  days  of  usable  data) 

Discharge  (ft3/s) 

-285 

96 

Nitrogen  Load  (short  tons) 

-12 

7 

-1.73 

-2.38 

Phosphorus  Load  (short  tons) 

-4 

3 

-1.34 

-1.83 

Sept.  7-Oct.  5,  1994 

Volume  (Mgal) 

186,133 

201,978 

.92 

1.26 

(24  days  of  usable  data) 

Discharge  (ft3/s) 

12,000 

13,022 

Nitrogen  Load  (short  tons) 

913 

1,276 

.72 

.98 

Phosphorus  Load  (short  tons) 

312 

539 

.58 

.79 

April  17-May  17,  1995 

Volume  (Mgal) 

-1,897 

12,433 

-  .15 

-  .21 

(26  days  of  usable  data) 

Discharge  (ft3/s) 

-113 

740 

Nitrogen  Load  (short  tons) 

-9 

91 

-  .10 

-  .14 

Phosphorus  Load  (short  tons) 

-3 

40 

-  .08 

-  .11 

Oct.  23-Nov.  21,  1995 

Volume  (Mgal) 

28,051 

33,612 

.83 

1.14 

(17  days  of  usable  data) 

Discharge  (ft3/s) 

2,553 

3,059 

Nitrogen  Load  (short  tons) 

138 

238 

.58 

.79 

Phosphorus  Load  (short  tons) 

47 

102 

.46 

.63 

Tor  1994,  Basin-ratio  estimation  =  (value  at  Station  1  +  value  at  Station  2  +  value  at  Station  3)  *  1.37. 

For  1995,  Basin-ratio  estimation  =  (value  at  Station  1  +  value  at  Station  2  +  value  at  Station  3  +  value  at  Station  4)  *  1.12. 
2For  1994,  Corrected  basin  ratio  *  (value  at  Station  1  +  value  at  Station  2  +  values  at  Station  3)  =  Tidal  estimation. 

For  1995,  Corrected  basin  ratio  *  (value  at  Station  1  +  value  at  Station  2  +  value  at  Station  3  +  values  at  Station  4) 

=  Tidal  estimation. 


mouth.  For  example,  in  figure  12b,  the  basin-ratio 
discharge  shows  a  storm  event  beginning  on  April  28 
and  peaking  on  April  30,  which  actually  is  when  the 
storm  was  measured  at  the  upstream  stations.  The  tidal- 
estimation  discharge  shows  the  storm  event  beginning 
on  April  30  and  peaking  on  May  2,  which  indicates  that 
there  was  about  a  2-day  travel  time  for  this  storm  peak 
to  move  down  the  channel  to  the  mouth  of  the  river. 


daily  mean  discharges  and  loads,  indicating  that  water 
and  nutrients  were  going  into  storage  in  the  reach  of 
the  Peace  River  upstream  from  the  mouth  and  down¬ 
stream  from  the  long-term  gaging  stations.  The  basin- 
ratio  method  will  always  produce  positive  daily  mean 
discharges  and  loads,  because  computations  are  based 
on  discharges  and  loads  at  the  nontidal-upstream 
stations  which  are  always  positive. 


Daily  mean  discharge  at  the  Peace  River  at 
Punta  Gorda,  computed  from  filtered  tidal-estimation 
data,  ranged  from  an  average  of  -285  ft3/s  during  May- 
June  1994  to  an  average  of  12,000  ft3/s  during  Septern- 
ber-October  1994  (table  4).  Daily  mean  discharge 
computed  using  the  basin-ratio  method  ranged  from 
96  ft3/s  during  May-June  1994  to  13,022  ft3/s  during 
September-October  1994.  During  the  1994  and  1995 
low-flow  data-collection  periods,  basin-ratio  daily 
mean  discharges  were  96  ft3/s  and  740  ft3/s,  respec¬ 
tively,  whereas  the  tidal-estimation  daily  mean  dis¬ 
charges  were  -285  ft3/s  and  -113  ft3/s,  respectively. 
During  the  two  low-flow  deployments,  computations 
using  the  tidal-estimation  method  produced  negative 


Total  discharge  and  total  nutrient  loads  com¬ 
puted  for  the  Peace  River  at  Punta  Gorda  for  each 
deployment  arc  presented  in  table  4.  Using  the  tidal- 
estimation  method,  total  nitrogen  loads  ranged  from 
-12  short  tons  during  May-June  1994  to  913  short  tons 
during  September-October  1994.  Using  the  basin-ratio 
method,  total  nitrogen  loads  ranged  from  7  short  tons 
to  1,276  short  tons  for  those  deployments,  respec¬ 
tively.  Using  the  tidal-estimation  method,  total  phos¬ 
phorus  loads  ranged  from  -4  short  tons  during  May- 
June  1994  to  312  short  tons  during  September-October 
1994.  Total  phosphorus  loads  using  the  basin-ratio 
method  ranged  from  3  to  539  short  tons  for  those 
deployments,  respectively. 
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Tidal-estimation  discharges  and  nutrient  loads 
were  less  than  basin-ratio  estimates  during  the  two 
high-flow  periods  of  September- October  1994  and 
October-November  1995,  which  corresponds  to  the 
results  for  the  Myakka  River  at  El  Jobean.  For  the  two 
high-flow  deployments,  tidal-estimation  discharges 
were  less  than  basin-ratio  discharges  by  factors  of  0.92 
and  0.83,  nitrogen  loads  were  less  than  the  basin-ratio 
nitrogen  loads  by  factors  of  0.72  and  0.58,  and  phos¬ 
phorus  loads  were  less  than  the  basin-ratio  phosphorus 
loads  by  factors  of  0.58  and  0.46.  The  differences 
between  the  tidal-estimation  and  basin-ratio  values  for 
nitrogen  and  phosphorus  loads  during  the  two  high- 
flow  periods  indicate  that  the  nitrogen  and  phosphorus 
concentrations  in  the  ungaged-tidal  basin  differ  from 
the  concentrations  in  the  gaged-nontidal  basin.  The 
ratios  of  tidal-estimation  method/basin-ratio  method 
for  discharges  were  always  greater  than  the  ratios  for 
nitrogen  and  phosphorus  load.  As  observed  for  the 
Myakka  River,  nutrient  concentrations  are  not  the 
same  throughout  the  Peace  River  Basin,  and  the  ratios 
of  tidal-estimation/basin  ratio  discharges  and  nutrient 
loads  are,  therefore,  not  the  same. 

For  the  two  low-flow  periods,  tidal-estimation 
discharges  and  nutrient  loads  at  the  Myakka  River  at 
El  Jobean  were  greater  than  basin-ratio  estimates.  In 
contrast,  at  the  Peace  River  at  Punta  Gorda,  not  only 
were  tidal-estimation  discharges  and  nutrient  loads 
less  than  basin-ratio  estimates,  the  daily  mean  dis¬ 
charges  for  the  low-flow  deployments  were  negative. 
Because  discharges  were  negative,  nutrient  loads  were 
also  negative,  indicating  that  water  and  nutrients  were 
going  into  storage  in  the  reach  upstream  from  the 
mouth  of  the  Peace  River  and  downstream  from  the 
long-term  gaging  stations. 

As  noted  in  the  discussion  for  the  Myakka 
River,  the  errors  associated  with  both  methods  tend  to 
increase  as  discharge  decreases.  The  daily  mean  dis¬ 
charge  during  the  low-flow  periods  (-285  to  -113  ft3/s) 
is  very  small  relative  to  the  magnitude  of  instanta¬ 
neous  dischai'gcs,  ±  40,000  ft3/s  and  to  the  average 
standard  error  in  measured  instantaneous  discharge, 

±  3,600  ft3/s.  The  measurement  section  near  the  mouth 
of  the  Peace  River  is  more  than  three  times  as  wide  as 
the  section  near  the  mouth  of  the  Myakka  River  and, 
consequently,  the  accuracy  of  the  index/mean  velocity 
rating  is  more  sensitive  to  the  location  of  the  index 
velocity  SIP.  Also,  discharge  through  700  ft  of  the 
4,100-ft  cross  section  was  estimated  because  of  the 
depth  limitations  of  the  ADCP. 


Potentiometric-surface  maps  of  the  intermediate 
and  Floridan  aquifer  systems  for  1994  and  1995  show 
a  potential  for  ground-water  discharge  into  the  tidal 
reach  of  the  Peace  River  (Metz  and  Brendle,  1994; 
Metz,  1995;  Metz  and  Stelman,  1995;  Metz  and  others 
1996).  It  is,  therefore,  improbable  that  surface  waters 
could  be  lost  through  recharge  or  go  into  storage  in  the 
ground-water  system. 

There  is  significant  potential  for  storage  in  the 
tidal  reach  of  the  Peace  River  upstream  from  U.S.  41. 
The  open-water  surface  area  of  the  tidal  reach 
upstream  from  U.S.  41  is  about  10  mi2  and  there  is 
substantial  additional  storage  in  adjacent  mangroves 
and  salt  marshes.  Considering  just  the  open- water  sur¬ 
face  area,  inflow  would  have  to  exceed  outflow  by 
32  ft3/s  over  the  course  of  1  day  for  stage  in  the  reach 
to  increase  0.01  ft.  Inflow  would  have  to  exceed  out¬ 
flow  by  320  ft3/s  over  the  course  of  1  day  for  the  stage 
in  the  reach  to  increase  0.1  ft.  The  13  days  of  data 
available  during  the  May-June  1994  deployment  rep¬ 
resent  a  tidal  range  that  began  during  a  neap-tide 
period  and  ended  during  a  spring-tide  period.  On  May 
18,  1994,  the  tidal  range  measured  at  Punta  Gorda  was 
1.5  ft  and  the  maximum  tidal  stage  was  1.35  ft  above 
sea  level.  By  the  time  of  the  full  moon  and  maximum 
spring  tide  on  May  25-26,  the  measured  tidal  range 
was  2.65  ft  and  the  maximum  tidal  stage  was  2.51  ft 
above  sea  level.  Thus,  the  hydrodynamic  conditions 
indicate  that  during  this  deployment,  the  amount  of 
water  moving  upstream  from  the  mouth  by  the  action 
of  the  tides  would  tend  to  be  greater  than  the  amount 
of  water  moving  downstream  and  would  cause  water 
to  go  into  storage  upstream  from  the  mouth. 

The  May-June  1994  deployment  on  the  Peace 
River  illustrates  several  problems  associated  with 
short-term  deployments.  The  available  data  from  the 
deployment  began  during  a  neap  tide  and  ended  during 
a  spring  tide.  If  the  available  continuous  data  spanned 
a  longer  period  of  time  it  would  be  possible  to  more 
accurately  identify  temporary  storage  effects.  The 
drawbacks  associated  with  using  the  basin-ratio 
method  for  short-term  estimates  of  discharge  and  load¬ 
ing  are  also  apparent  because  the  basin-ratio  method 
cannot  account  for  the  time  lags  associated  with  vari¬ 
ous  tributaries  or  for  the  temporary  storage  effects  that 
are  associated  with  the  transition  from  one  tidal  range 
to  another. 
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SUMMARY  AND  CONCLUSIONS 

Nutrient  enrichment  has  been  identified  as  a 
significant  concern  in  the  Charlotte  Harbor  estuarine 
system,  and  consequently,  it  is  important  to  accurately 
estimate  the  magnitude  of  nutrient  loads  transported  by 
tributaries  to  the  estuary.  The  Myakka  and  Peace  Rivers 
constitute  more  than  60  percent  of  the  total  inflow  area 
and  contribute  more  than  half  of  the  total  inflow  to 
Charlotte  Harbor. 

Two  methods  for  estimating  discharge  and  nutri¬ 
ent  loads  from  tidally  affected  reaches  of  the  Myakka 
and  Peace  Rivers  were  compared.  The  first  method  was 
a  tidal-estimation  method,  in  which  nutrient  loads  were 
estimated  based  on  stage,  water- velocity,  discharge,  and 
water-quality  data  collected  near-  the  mouths  of  the  two 
rivers.  The  second  method  was  a  traditional  basin-ratio 
method  in  which  discharge  and  nutrient  loads  at  the 
mouths  were  estimated  from  discharge  and  loads  mea¬ 
sured  at  upstream  stations.  Data  were  collected  during 
four  independent  30-day  periods  during  1994  and  1995. 
The  measured  data  were  filtered  to  remove  short-term 
(less  than  2  days)  variations  resulting  from  mixed  tides, 
and  the  filtered  data  were  used  to  estimate  daily  mean 
discharges,  nitrogen  loads,  and  phosphorus  loads  for  the 
two  rivers.  Values  computed  using  the  tidal-estimation 
method  were  then  compared  to  values  computed  using 
the  basin-ratio  method. 

Several  factors  affect  the  relative  accuracy  of 
each  of  the  methods  and  must  be  considered  in  deter¬ 
mining  which  method  produces  estimates  that  are  most 
representative  of  actual  discharges  and  nutrient  loads. 
The  basin-ratio  method  assumes  that  the  quality  and 
quantity  of  runoff  and  base  flow  are  the  same  in 
ungaged,  tidal  parts  of  the  basin  as  they  are  in  the 
gaged,  nontidal  part  of  the  basin.  However,  topography 
and  land  use  are  substantially  different  in  the  ungaged, 
coastal  parts  of  the  Myakka  and  Peace  River  Basins 
than  in  the  gaged,  upland  parts.  Also,  the  basin-ratio 
method  presumes  that  discharge  is  hydraulically  con¬ 
trolled  at  the  upstream  end  of  the  reach.  In  the  Myakka 
and  Peace  Rivers  water-surface  elevations  at  the  down¬ 
stream  end  of  the  tidal  reaches,  which  respond  to  astro¬ 
nomical  and  meteorological  forces,  significantly  affect 
discharge.  The  tidal-estimation  method  is  subject  to 
measurement  errors  that  can  be  significant  relative  to 
the  magnitude  of  daily  mean  discharge. 

Estimates  of  discharge  and  nutrient  loads  com¬ 
puted  from  the  tidal-estimation  and  basin-ratio  methods 
were  most  similar  during  high-flow  periods.  However, 
during  high  flow,  the  values  computed  from  the  tidal- 


estimation  method  were  consistently  lower  than  the 
values  computed  from  the  basin-ratio  method.  The 
Myakka  River  discharges  computed  from  the  tidal-esti¬ 
mation  method  were  lower  than  those  from  the  basin- 
ratio  method  by  factors  of  0.77  and  0.74.  The  Peace 
River  discharges  computed  from  the  tidal-estimation 
method  were  lower  by  factors  of  0.92  and  0.83.  Nitro¬ 
gen  loads  computed  from  the  tidal-estimation  method 
were  lower  than  those  from  the  basin-ratio  method  by 
factors  of  0.74  and  0.57  for  the  Myakka  River  and  by 
factors  of  0.72  and  0.58  for  the  Peace  River.  Phosphorus 
loads  were  lower  by  factors  of  0.60  and  0.33  for  the 
Myakka  River  and  by  0.58  and  0.46  for  the  Peace  River. 

During  low-flow  periods,  there  were  substantial 
differences  between  discharges  and  nutrient  loads  com¬ 
puted  from  the  tidal-estimation  and  basin-ratio  meth¬ 
ods.  Furthermore,  the  differences  between  the  methods 
were  not  consistent.  The  Myakka  River  discharges 
computed  from  the  tidal-estimation  method  were  higher 
than  those  computed  from  the  basin-ratio  method  by 
factors  of  2.88  and  5.22.  Nitrogen  loads  were  higher  by 
factors  of  2.18  and  3.89  and  phosphorus  loads  were 
higher  by  factors  of  1.31  and  2.34.  The  Peace  River 
discharges  and  nutrient  loads  computed  by  the  tidal- 
estimation  method  were  not  only  lower  than  those  com¬ 
puted  from  the  basin-ratio  method,  but  they  actually 
reflected  a  negative,  or  upstream,  net  movement. 

During  the  1994  low-flow  data-collection  period,  daily 
mean  discharge  at  the  mouth  of  the  Peace  River,  com¬ 
puted  by  the  traditional  basin-ratio  method  was  96  ft3/s. 
Using  the  tidal-estimation  method,  daily  mean  dis¬ 
charge  was  -285  ft3/s.  Hydrodynamic  conditions  for  the 
13  days  of  available  data  during  this  deployment  indi¬ 
cate  that  the  amount  of  water  moving  upstream  as  a 
result  of  tidal  forces  would  tend  to  be  greater  than  the 
amount  of  water  moving  downstream.  During  the  1995 
low-flow  data-collection  period,  basin-ratio  daily  mean 
discharge  was  740  ft3/ s  and  tidal-estimation  daily  mean 
discharge  was  -113  ft3/s.  Because  nutrient  loads  are  a 
function  of  discharge,  the  basin-ratio  method  produced 
positive  net  loads  and  the  tidal-estimation  method 
produced  negative  net  loads. 

Negative  net  loads  could  be,  in  part,  the  result  of 
measurement  error  as  well  as  representing  water  and 
nutrients  going  into  storage.  Short-term  tidal  measure¬ 
ment  results  should  be  used  with  caution,  because 
antecedent  conditions  can  influence  the  discharge  and 
nutrient  loads.  Continuous  tidal  data  collected  over  a  1- 
or  2-year  period  would  be  necessary  to  more  accurately 
estimate  the  tidally  affected  discharge  and  nutrient  loads 
for  the  Myakka  and  Peace  River  Basins. 
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